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■ Abstract: We compute for various perturbed conformal field theories the vacuum en- 

ergies by means of the thermodynamic Bethe ansatz. Depending on the infrared and 
ultraviolet divergencies of the models, governed by the scaling dimensions of the under- 
lying perturbed conformal field theory in the ultraviolet, the vacuum energies exhibit 



o 
^1- . 

different types of characteristics. In particular, for the homogeneous sine-Gordon models 
we observe that once the conformal dimension of the perturbing scalar field is smaller 
or greater than 1/2, the vacuum energies are positive or negative, respectively. This be- 
haviour indicates the need for additional ultraviolet counterterms in the latter case. At 
• the transition points we obtain an infinite vacuum energy, which is partly explainable 



with the presence of several free Fermions in the models studied. 



1. Introduction 

According to the ideas developed first in Q a large class of massive quantum field theories 
in 1+1 space-time dimensions can be viewed as perturbed conformal field theories (CFT) 
with Euclidian action 

5 = 5 C FT + A d 2 xip(x) . (1.1) 



Here 5cft denotes a fixed point action, (p(x) a scalar field with (left, right) conformal 
dimension (A, A) and A a coupling constant, scaling with (1 — A, 1 — A). The great virtue 
of such theories is that very often they are integrable and can be solved exactly, that is to 
all orders in perturbation theory. Since the original formulation various non-perturbative 
techniques have been developed to study such theories with great success. Nonetheless, 
once the CFT is well investigated one may also employ standard perturbative arguments 
and unravel the meaning of certain types of behaviour in that more traditional language. 
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Accordingly, the vacuum expectation value of any local operator O can then be computed 
as 

(0(z,z)) = Z^J^^-T / d 2 z 1 ...d 2 z n (0(z,z)<p(z 1 ,z 1 )...ip(z n ,z n )} CFT . (1.2) 

n=0 ' J 

Here the normalization factor is in general Z = <^exp — A J d 2 z(p(z, z)} CFT , with ) CFT de- 
noting the vacuum state related to 5cft- In quantum field theories such expressions are 
plagued by various types of divergencies. First of all one encounters the infinities due to 
the self-contraction of the fields, which can be regularized fairly easily by a normal ordering 
prescription. Second, one might have ultraviolet (UV) singularities for (z — z^) — > 0. Here 
the case O =ip will be important, for which we can approximate with the help of standard 
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CFT operator product expansion the integrals in (1.2) as ~ f dzi \ z — Zi\ . Thus for 



A < 1/2 the integrals in (1.2) remain finite, whereas for A > 1/2 we require in general 
counterterms to elliminate the divergencies. Third, one might have infrared (IR) singular- 
ties for (z — Zi) — > 00. In the infinite plane it is usually an intricate issue to handle them 
B ||]. However, when formulating the theory from the very beginning on a cylinder instead 
of an infinite plane the integrals in ( |1.2| ) will automatically be IR finite for A > 0, as the 
cylinder radius R constitutes a natural cut off. The fourth singularity occurring is related 
to the fact, that even when the individual integrals in (|l.2j ) are finite the entire series will 
in general be IR divergent for large R. 

Supposing now that one is able to compute ([h^) exactly, that is to all orders in 
perturbation theory, the different types of renormalization quantities should be tractable 
in that context. In fact, the thermodynamic Bethe ansatz (TBA) Q is a method which 
allows such identifications when O is taken to be the energy operator. The above mentioned 
arguments hold when recalling || that this operator is proportional to the perturbing field 
ip. Defining then for the ground state energy E(R) the scaling function c(R) = —6RE(R)/ir 
one encounters several types of general behaviours, which can all be brought into the generic 
form 

00 

c(r) =c eS + £ r 2 +£' r 2 lnr + Y. £ " Xn f^ ■ ( L3 ) 

n=l 

Usually one uses the dimensionless parameter r = R/m with m being a mass scale and 
£0, £q being finite real numbers. The function c(r) is normalized in such a way that 
c(r = 0) coincides with the effective central charge c e ff = c — 24A m i n , with c being the 
Virasoro central charge of the underlying ultraviolet conformal field theory and A m j n the 
smallest conformal scaling dimension in the model. This constant c has the well known 
interpretation as the Casimir energy, which is the vacuum energy on the cylinder and 
becomes zero when mapped onto the plane. Viewing (|1.2[) as resulting from a partition 



function, the term £qt has to be present in (1.3), since thermodynamics dictates that 



for large r the energy has to be proportional to the volume. In quantum field theoretic 
terms both £qt 2 and £^Qr 2 lnr are related to renormalization issues, characterized by the 
conformal dimension A as described above. These terms are also needed in order to ensure 
that linv^oo c(r) = 0, which one expects for a purely massive model. Finally, the f n (f) 
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result from the integrals in ( |1.2[) and takes on various general forms depending on the 
regime in which A is valued. 

In this paper we will discuss more concretely the precise nature of the expansion (|1.3|). 
We will first recall in section 2 how the TBA can be used to compute the vacuum energies 
and in the following sections we discuss the different regimes for different types of concrete 
theories, the homogeneous sine-Gordon (HSG) models || [?]] and affine Toda field theories 
(ATFT) [|. These theories probe several regimes for A and exhibit different types 
of behaviours. In particular for the HSG- models, which are defined in the entire regime 
< A < 1, our results will be new. Our conclusions are stated in section 6. 



2. Vacuum energies from the TBA 

Let us briefly recall the main steps of how vacuum energies may be computed Q (more 
details on the arguments can also be found in [pX|]) non-perturbatively with the help of the 
TBA. One considers a relativistic theory in which the scattering matrices Sij{6) for the 
particles of the type i,j with masses mi,nij are known as functions of the rapidity 9. Then 
the entire TBA analysis can be formulated with only two inputs: first the dynamical inter- 
action, which enters via the logarithmic derivative of the S-matrix <fij(9) = —idlnSij(9)/d9 
and an assumption on the statistical interaction, which we choose here to be of fermionic 
type. The thermodynamic Bethe ansatz equations are then a set of coupled non-linear 
integral equations 



cosh 9 = Ei(e, r) + Y^ifij * ln(l + e -e ')](0, r) , (TBA) 



where the pseudo-energies 6i(9,r) are the unknown quantities. We denote as usual the 
convolution of two functions by (/ * g) (9) := 1/(2-7!" ) J d9'f(9 — 9')g{9'). The scaling pa- 
rameter is related to the inverse temperature T as r = m/T, mi — > mi/m, with m being 
an overall mass scale. In Q it was shown that when taking the sum and difference of the 
derivatives d/dr ^TBA ) and d/d9 ( [TBA )/r one may derive a set of coupled linear integral 



equations for the quantities 



A( ,, r) _*^) ± i*&l2, (2 ,) 



respectively, namely 

^4(0, r) = m t e ±e + J> y * -J—^ ± ](9,r) . (2.2) 



The strategic is now to solve first the equations ( [TBA| ) for Ei(9,r) and thereafter (2.2) for 



ip±{9, r). Once one has carried out the first step, one can already compute the scaling 
function 

oo 

c(r) = ^j^2 m i J d9 cosh 9 L { (9, r) , (2.3) 
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with Li(0,r) = ln(l + e~ £i ( 9 ' r ^). Concerning the status of analytical solutions for (|2.2|), it is 
similar as for the TBA itself, that is only for free theories a closed solution was found 
and for interacting theories was only solved in the extreme ultraviolet limit. Numerical 
solutions exist even less. Once it is solved, one may compute the vacuum expectation value 
of the trace of the energy momentum tensor, i.e. vacuum energies 



i 

oo 



,r e 



(2.4) 
(2.5) 



In a parity invariant theory we have Ei(0,r) = e$(— 0, r) and consequently ifj\(9,r) = 
ip l _(9,r), TI = T % _ = T l such that matters simplify. We like to keep the treatment here 
generic for a while as we will also consider below the homogeneous sine-Gordon models, 
which are not parity invariant. 

There exists no systematic way to solve the equations ( [TBA ) and ( |2.2| ) analytically, 
albeit, numerically this is a solvable problem. Nonetheless, it is well known that at the fixed 
points approximations can be made, such that one can solve ( [TBA[ ) analytically and hence 
also obtain analytic expressions for (|2.3| ) at these points (r = is one of them) . Likewise we 



expect to be able to solve (2^2) for these values and compute (T^) analytically. Following 



now essentially the argumentation of [||, |T^], we need to make only three assumptions: 



i) The logarithmic derivative of the scattering matrix in (TBA) admits an expansion of 
the form 

^) = -J2^ e ~ sW - ( 2 - 6 ) 



ii) For the first coefficient in ( |2.6| ) we presume proportionality to the masses 

Vij = Pij m i m j 
for some function specific to the particular theory. 

iii) One assumes that 



§i(0) -e< < e 6 



for < 



(2.7) 



(2.8) 



where the £j are the pseudo-energies of the constant TBA equation and the £i(9) are 
quantities in the r-independent TBA-equation 



(fij * Lj(0) = -EiiO) + mie 



(2- 



obtained from ( [TBA ) by the shift 6 — * + ln(r/2), £i(6,r) — > £j(0). This assumption 
is usually difficult to justify a priori, but is sustained in hindsight by meaningful 
results or supported by numerical data. 
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For 6 —* — oo one can now derive with (|2.8| ) the equation 



Vij * L 3 (9) = -e % + ^VgM + 0(e" e ) (2.10) 



where r\ > 2. Comparing then ( |2.9| ) and ( 2.10| ) for the parity invariant case, it follows 



directly with (2.8) that 

m . = — J$Ti . (2.11) 

* 27T %3 

Finally we deduce the expression for the vacuum expectation value for the energy momen- 
tum tensor with ( |2~8| ) and ([O]) to 

<n> = ^E^- 1 - ( 2 - 12 ) 

This quantity is of course sensitive to above mentioned renormalization issues and possibly 
exhibits the distinction between the different regimes quoted. Furthermore, one has the 
possibility of comparison, as there are various other methods to obtain the vacuum energies, 



such as the truncated conformal space approach [11] or a matching between the high-energy 
behaviour of the scattering matrix with a Feynman diagramatic analysis |l2j . 

Let us briefly comment on the different regimes: 
< A < 1/2 : As mentioned in the introduction, in this regime the individual integrals 
in the expansion ([l]^) are UV and IR convergent term by term when formulated on the 
cylinder. From general arguments one finds for the behaviour in ( |1.3[) that £ = and 



f n (r) = r 2n(i-A) g_ From and (gj) follows also that we can identify (T^) \ r=0 = 
—tt 2 /2>£q. Thermodynamically this term can be seen as the infinite volume energy and field 
theoretically this corresponds to the sum of all infrared substractions, which achieve the 
convergence of the sums ([L.3|) for large r. 



1/2 < A < 1 : Now the individual integrals in the expansion Ql.2| ) are still IR convergent, 
but cease to be UV convergent. Nonetheless, we may still employ similar arguments as 
in the previous regime and find again for the behaviour in ( |1.3| ) that £' = and f n ( r ) = 
r 2n(i-A) From (|1.3|) and (2.4) follows once more that we can identify (T^) | r =o = 



' [IJ-OIJ. 1 iUlli U-L.OI/ ailU. U^.'ii/ lUHUWO UllV^C IIIUIC blLfXb WC lU-Clltll^ \J- f+ 

—tt 2 /3£(). However, now the field theoretic interpretation of this term changes. Since 
we require in this case UV conterterms to make the individual integrals finite, the £q- 
term corresponds now to the sum of these UV counterterms and all infrared substractions, 
which achieve the convergence of the sums ( |1.3| ) for large r. Indeed, for the concrete models 
studied below this becomes visible in a change of sign in the transition from the regime 
A < 1/2 to A > 1/2. 

A < : Now the individual integrals in the expansion ( |1.3[ ) are still UV convergent, but 
cease to be IR convergent even on a cylinder. General arguments now yield for the be- 
haviour in ( |1.3| ) that £q ^ 0, £' Q ^ and / n (r) = (a + \n(r))~ n with a being some constant 
ig, [L§ 0, 0, One still finds that (T^) | r=0 = -tt 2 /3£ , e.g. 0], but now the 

interpretation is less obvious as some counterterms also accumulate in the £ -term. 
A = 1/2 : In this case one usually finds free Fermions in the model and £q ^ 0, £' ^ 0, 



fn( r ) = T n . Now the vacuuum energy is divergent, see e.g. [1C] for an analytical expression. 
We will investigate some concrete theories. 
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3. < A < 1/2, minimal affine Toda field theories 

These theories have been studied before flO| , |l9|| , nonetheless, we recall them here as they 
are easy examples which illustrate the working of the above formulae and we shall also 
point out some novel features. We recall first that minimal affine Toda field theories can 
be realized as perturbations of the coset conformal field theories gi®gi/g2> with g& being a 
simply laced Kac-Moody algebra of rank i and level k [20, ETA. The corresponding Virasoro 
central charges c and conformal dimension of the perturbing operator A are 

21 2 

and A = , (3.1) 



2 + h 2 + h 

respectively. Apart from h = 2, i.e. the free Fermion with g = Ax, we always have 
for the Coxeter number h > 2 and are therefore in the stated regime < A < 1/2. 
The renormalization issues are handled most easily in this case and the vacuum energies 
are computable with the above arguments. With regard to assumption i), we recall the 
expansion of the TBA-kernel for these theories |l0|, ^2|, [23j 

Vij (9) = -i^cot^x^x^e-^ , (3.2) 

with £ ={s + nh}, s being an exponent of g, n G No and Xi(s) are the left eigenvectors of 
the Cartan matrix. In particular, we have Xj(l) = rrii/m, with m being an overall mass 
scale, which is needed for the assumption ii) to hold. Having therefore the quantity <p\j in 



the form ( |2.7| ), we deduce immediately with the help of (2.12) 



7T 7T 



(T^) = m 2 - tan - . (3.3) 

Obviously apart from the free Fermion with h = 2, when (T^) — > oo, we have (T M M ) > 0. 
This result agrees with a similar formula obtained in in terms of the coefficients ip^ 
without explicit evaluation and "1" referring to the lightest particle. More concret case- 
by-case studies were carried out in [T^] for perturbations of g/ (g> gfc/gfc + ;-coset CFT's (see 
formulae (3.14) therein). When using the overall mass scale to perform suitable normal- 
izations the formula for k = I = 1 in there can be brought into the universal formula (|3.3;), 
which is not obvious at first sight. The formulae in fT^j are expressed in terms of a mass 
scale M whose relation with our m varies for every theory as 

At : M = m sin 
D e : M = m/y/2 



E 6 : M = m\J § sin § (3.4) 



E 7 : M = my'sin^/sin^ 
E 8 : M = my^sin^jsinf . 

The advantage of our formulation relies on the fact that the masses are normalized with 
respect to the same general mass scale m for all simply laced Lie algebras, which allows 
for the very compact and generic expression (^^). Alternatively these results were also 
confirmed in [24]. 
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4. < A < 1, g fc -homogeneous sine-Gordon models 

Let us now consider a theory which is more interesting with regard to the above mentioned 
problematic, namely the g^-HSG model |(], M, with g being a simple Lie algebra of rank 
£ and level k. These models can be viewed as perturbed Wess-Zumino-Novikov-Witten 
(WZNW) §§ coset-models 

Srsg = Swznw + 52 / ^ x ( A +>9(x)~ A-g(x)) . (4.1) 
up J 

Here ( , ) denotes the Killing form of g and g(x) a group valued bosonic scalar field. A± 
are arbitrary semi-simple elements of the Cartan subalgebra associated with the maximal 
abelian torus h C g, which have to be chosen not to be orthogonal to any root of g. The 
parameters m and (5 are the bare mass scale and the coupling constant, respectively. The 
Virasoro central charge of the coset model and the dimension of the perturbing operator 
are computed to 

n kh-h y , /i v 

c = £ — and A = — , (4.2) 

k + h v k + h y K ' 

with (/i v ) h being the (dual) Coxeter number of g. We note that now the constraint A < 1/2 

does not automatically hold for each level and the above mentioned complications could 

arise for some theories in this series when changing from k > h v to k < h v . Up to now no 

indication for a different behaviour of the theories in this two different regimes have been 

found in the literature. We treat the simply laced and non-simply laced cases separately. 

4.1 Simply laced HSG-models 

As in the original formulation of these models, the algebra g is assumed to be simply laced. 
Since for this case the expansion of the kernel <p does not appear in the literature, we will 



start with the scattering matrix, which was found originally in []2J|] (see [27] for an integral 
representation). We cast the matrix describing the scattering between the particle of type 
A = (a, a) and B = (6, b), with 1 < a,b < £; 1 < a,b < k into the form 

s%«» = d » P / f^ f^^Lt m ^"' ■ < 4 ' 3 > 

Here r]^ = exp[i7re-j ) (2 — lA k _ 1 )^, \ are constant phase factors not relevant for our analysis, 
^Sfc(^) = 25-^ cosh t/A; — 1-^ with / being the incidence matrix of g and the cr's are the 
resonance parameters, which indicate the presence of unstable particles in these models. 
In order to evaluate the expansion for ip, we can employ the residue theorem for a contour 
along the real axis closing up in the positive half of the complex plane encircling all poles 
on the imaginary axis in the upper half plane. Noting that in fl4.3| ) t = inn are simple 
poles, except for t = innk which constitute double poles for n G N, we deduce for a-^ = 

4(A)=W V Res k {t) ^m^zMK e ^ (4.4) 

^ abK 1 ^ t=i7Ts\ it J ab J sinh(t/A;)sinht v ' 

s=l;s^nk V 7 V ' 1 

= _ 2 £ 4 iW "^yt) e -.i,i . (4 . 5) 

s=l;s^nk ^ ' ' 
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The desired coefficient ( |2.7D follows from this directly to 



(i) 



^ab 



2K^{in)m a mb/m 2 sin(-7r/£;) 



(4.6) 



where = m a m a with m a = sina7r//c being the masses of Afc_i-affme Toda field theory 
and m a are I free mass scales. We choose them here to be all equal m a = m V a. Finally 
we derive from this a closed expression for the vacuum expectation value for the trace of 
energy-momentum tensor 



Km 



sin(7r/A;) 



ab 



(4.7) 



a,6=l 



We are not aware of a generic formulation for K 1 (ik) and analyse therefore the expression 
Q4.7|) below in more detail case-by-case. We can summarize our findings as 



> for k > h = A < 1/2 
-> oo for k = h = A = 1/2 
< for k < h = A > 1/2 



(4.8) 



In many cases we can attribute the divergence for A = 1/2 to the presence of free Fer- 
minons. The change of sign when going from A<l/2toA>l/2 reflects the fact that 
besides the IR counterterms, which achieve the convergence of the sums ( |1.2|) for large r, 
needed in both cases in the latter we also require UV conterterms to make the individual 
integrals in (|1.2|) finite. 



4.1.1 (A £ ) fc -HSG model 



For An the Coxeter number is h = £ + 1. The inverse of the matrix relevant in (4.7) can be 
cast in this case into a simple formula 



sin(a7r/fc) sin[(/i — b)ir/k] 



for a < b. 



ab s'm(ir/k) sm(hir/k) 

Computing the sums over both entries then yields after some algebra 



-Km 



2 tan 2 K/2k 
Hence, the condition (T^) > becomes 



YlK 7T 

tan h tan — 

2k 2k 



Kk 7T 

tan — - > h tan — 
2k 2k 



or equivalently, when expanding the tan, 

1 



4/i^ 

7T k. 



,41 

>h--y ■ 

TT k ^-^ 



{l/kf 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



vr k ^ (2n - l) 2 - (h/k) 2 vrfc^(2n-l) 2 

It is easily seen that ( |4 . 1 2| ) holds term by term once h/k < 1, hence establishing the first 
inequality in (|4.8| ) . Similar arguments show that the opposite inequality holds in the regime 
h/k > 1. We comment more on the case k = h below. 
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4.1.2 (D^fc-HSG model 

For the Coxeter number is h = 2i — 2 and by evaluating (fOj) similarly as in the previous 
subsection, we find 



vrm 2 sin £ [2 - (2 + h) cos |f ] + 2-Kin 2 sin 



hm 
2k 



Sm 2k COS 2k 



<n> = fc 1 sin2 ; - z 1 ' ' ■ ' "* ' 

The condition (T^) > is now equivalent to 



7T 

sin — 

k 



(2 + h) 



<2tan— . (4.14) 
2k 



Expanding the left and right hand side of this inequality yields by similar arguments as in 
the previous subsection once more the relation (f4.8j). 

4.1.3 (E 6 ) fc -HSG model 

For Eg the Coxeter number is h = 12 and we find 



We see that the numerator is < for k = 2 and > for k > 2. The denominator is > for 
A; = 2, A; > 12 and < for 2 < k < 12. The denominator vanishes for fc = 12. Hence the 
relation ( fQj| ) holds. 

4.1.4 (E 7 ) fc -HSG model 

For E7 the Coxeter number is h = 18 and we find 
irm 2 y^n-i T v smpir/k 

T m = /ta « 71 9^ f = ( 9 > 18 > 20 > 22 ' 17 ' 12 > 7 ) • ( 416 ) 

p ' COS 7T//c(4 cos 6717 K — 2) 

We observe now that the numerator is < for k < 4 and > otherwise. The denominator 
on the other hand is > for k = 3, k > 18 and < otherwise except for k = 2, 18 in which 
case it is zero. Hence (|4.8|) holds also in this case. 



4.1.5 (Eg)fc-HSG model 

For Eg the Coxeter number is h = 30 and we find 
irm 2 y\ 1 t„ sm.pi: Ik 

= o 77- I ," ' 77 T70 ^ = ^ 8 ' 12 ' 12 ' 13 ' 10 ' 7 ' 4 ) • ( 41? ) 

p ' cos 87T/K + cos bir/k — cos 2ir/k — 1/2 

We see that the numerator is < for k < 5 and > otherwise. The denominator is > 
for k = 2, 3, 4; A; > 30 and < otherwise except for k = 30 in which case it is zero. Hence 
Q4.8D holds also in this case. 
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4.2 Non-simply laced HSG-models 

Now we allow the algebra g to be also non-simply laced. In this case the scattering matrix 
is slightly more complicated as the symmetry between the exchange of long and short roots 
is lost. It can be restored by the use of the symmetrizers t a of the incidence matrix of g, 
i.e. ta,I~r = trlr-, with t a = 2/a| and the length of long roots normalized to af = 2. In [28] 
the scattering between the particle of type A = (a, a) and B = (b,b), with 1 < a, b < £; 
1 < a, b < k a = t a k was proposed to be described by 

cab(^ _ ^ab pvn [ dt~ sinh(q£/fc a ) sinh[(l - b/k- b )t] it(e+ )/w 
S ab (e)-n ab expJ jK-Jt) sinh(tAa . )sinht e . (4.18) 

Here rf a h h are once more constant phase factors not relevant for our analysis. Further- 
more, one needs the quantity k ah - = femax(tj^j) and the matrix K with entries K~r(t) = 
25-^ cosh t jk a — I a }t~ h l max(ta j t^). Note that in comparison with p3] we interchanged the 
long and short roots, that is we have taken the fs to be left and not the right symmetrizers 
of the incidence matrix. A similar analysis as in the previous subsection leads now to the 
following expansion of the TBA-kernel 



fab 

such that 



SO = t ^(^) -(-/^^/^) e -^l, (4,9) 



s=l;s?nk ai M^/k- bJ 

tpV = 2K ai {iK)mlm h b /m 2 sm(n/k &i ) . (4.20) 

Here the masses are assumed to renormalize with an overall factor and are therefore ex- 
pected to be the same as in the semi-classical analysis ]29| ], that is m a a = msina-7r/k a . The 
overall mass scales associated with each colour are once more choosen to be the same. Thus 
we finally deduce 

I 

<r^> = ™ 2 £ Ik- 1 ] . (4.2i) 

— L A ab 

a,b=l 

where = K^i-n) sin(7r/A;-^). As in the previous case, we are not aware of a generic 
formulation for K' 1 and analyse therefore ( 4.21] ) in more detail case-by-case. Our findings 



are summarized as 

' > for k> h y ee A < 1/2 
^oo for k = h v = A = 1/2 , (4.22) 

< for k < h v = A < 1/2 

with similar interpretations as in (4.8). We establish ( |4.22j ) in more detail case-by-case. 



4.2.1 (G 2 ) fc -HSG model 

The dual Coxeter number for G2 is h v = 4 and the symmetrizers are taken to be t\ = 
3, £2 = 1- With these data we compute from (|4.2l|) 

, „ . 9 sin7r/A; + sin47r/3A; 

(T fl ..) = 2vrm 2 '— . , . (4.23) 

x M/ 2cos4vr/3/c - 1 v ; 
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Obviously, the numerator is > for k > 2, whereas the denominator is < for k = 2, 3 
and otherwise > except for k = 4 when it is zero. Evidently this agrees with ( gjg) . 



4.2.2 (F 4 ) fc -HSG model 

The dual Coxeter number for F4 is h v = 9 and the symmetrizers are taken to be t\ = ti = 1 
and t% = £4 = 2. From ( 4.21| ) we compute 



= 2E^ lB inW2fc-Bin37r/2fc 
v m/ 2cos3tt/A;-1 



(4.24) 



The numerator is > for k > 2, whereas the denominator is < for 2 < < 9 and 
otherwise > except for k = 9 when it is zero. Evidently this agrees with ( gUp . 



4.2.3 (B^) fc -HSG model 

The dual Coxeter number for is /i v = 2^ — 1 and the symmetrizers are taken to be 
ti = t% = . . . = tg-i = 2 and t# = 1. We find now for even rank £ 



h v — 1 (•£—■ 2)/2 

E sinf + fsin^ + 2cosf E (£ - 2p - 1) sin ^±£^1 



2 — 2fc ^ 2 0111 2fc 

7rm p=2 



p=i 



cos 



2k 



l + 2Ep / =i(-l) p cos. 



(4.25) 



whereas for odd I we obtain 



h v — 1 (£—3)/2 

, E sinf + fsin^ + 2cost E (/ - 2p - 1) an 
7rm p=i 



X 1 fi/ 



P =l 



cos 



2 k 



l + 2Ep / =i(-l) p cosf 



(4.26) 



Once more we confirm (|4.22j ). As the details are rather involvolved we drop them here. 



4.2.4 (Q) fe -HSG model 

The dual Coxeter number for Ci is h y =1 + 1 and the symmetrizers are taken to be 
t\ = t2 = ■ ■ ■ = = 1 and tt = 2. 



/rpfj, \ 



/rpfX \ 



irm 2 (i) 



cos 



2k 
2 



COS 



2k 



' E^b-^^l + Isin^ 
. l + 2E^i(-l) p cosf 

(p- l)sm — + -sm 

p=2 



2k 2 2k 



for t even 



for I odd. 



(4.27) 
(4.28) 



Once more we confirm ( 4.22 ) and drop the details for the same reasons as in the previous 
subsection. 
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5. A = 1/2, g h w -homogeneous sine-Gordon model 

The case A = 1/2 is very special as then the vacuum energy diverges in the extreme 
UV limit. Such type of behviour is well known from free Fermions in form of logarithmic 
ultraviolet singularities, meaning that fl2.4| ) yields (T^) — > oo for r — > 0. Explicit analytic 
formulae for the free Fermion c(r)-function can be found in pOj ]. Indeed in many cases we 
can make this connection quite explicit. It suffices to present an examples to illustrate this 
point. 

5.1 (^) £+1 -HSG theories 

Let us have a closer look at the (^Lj)£ + i-HSG theories in order to see how the Fermions arise 
in there. Obviously for h = k the expression ( EL~10l) yields (T^) -> oo. Already in @ [t 

was 



noticed that the (^42)3-HSG model decomposes into four free Fermions when the resonance 
parameter vanishes. From the fact that the central charge (^2) becomes in general £ 2 /2 



for (^)^_|_i-HSG models, one might suspect that they always decompose completely into 
£ 2 free Fermions for vanishing resonance parameters, such that each Fermion contributes 
1/2 to the central charge. However, this is not quite the the following argument 

shows. 

In order to count the Fermions, identified here simply with the amount of particles 
which contribute 1/2 to the central charge, we recall the constant TBA equations, which 
arise from ( [TBA ) after some standard manipulations. For the (A^+i-HSG models they 



take on the form 



x a 



J] (1 + 4)^ with N% = 5 ab 5- ai - (K Ae ) ab . (5.1 



Solving these equations for the x° =exp(— e") yields the effective central charge as 

— ££<i£iH (5 ' 2) 

a,a=l v a/ 

with C(x) = Y^n=x x n /n 2 + lnxln(l — x)/2 denoting Rogers dilogarithm. The solutions of 
( |5.1| ) are very simple in this case 

s _ sin[7ra/(l+l)] 
a ~ sin[7ra/(l+^)] • 

Therefore we have = x^ +1 ~ a = 1 and since £(1/2) = 7r 2 /12 it follows from this that 
each of the particles (a,a),(a,£ + 1 — a) for 1 < a < £ contributes 1/2 to the effective 



central charge in (5.2). Hence in the (j4.£)^+i-HSG models we have always 2£ or 2£ — 1 free 
Fermions when £ is odd or even, respectively. The remaining particles can be organized 
without exceptions in pairs (a, a), (a, a). Noting with ([5l]) that obviously x a a = (x?) _1 and 
recalling the fact that C(x) + C(l — x) = ir 2 /6 explains then that the central charge has to 
be an integer or a semi-integer for these models. 

In general, it is less straightforward for the other algebras to identify particles which 
directly contribute 1/2 to the central charge. In fact, mostly the particles occur in pairs, 
triplets or higher multiplets contributing integers or semi-integer values to c. 
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6. A < 0, affine Toda field theories 

Affine Toda field theories related to simply laced and non-simply laced Lie algebras have a 
quite different behaviour due to the fact that in the first case all masses renormalize with 



an overall factor, which is not the case in the latter (see e.g. [30]). As a result of this, the 
strong-weak duality observed for ATFT related to simply laced algebras is broken for those 
associated with non-simply laced Lie algebras. Despite the fact that there exists a uniform 
formulation, we will treat them here separately as this will be more transparent. 

6.1 Simply laced 

ATFT are quite well studied examples of integrable models, which can be viewed in the 
spirit of ( |l.l] ) which was noted first in [21, 20 1 



Satft= [ d 2 xhd^) 2 + ^m^ Si ' ■ (6-1) 



i=0 

The fixed point part of the action Soft corresponds to the conformal Toda field theories 
when the sum over the simple roots on starts at i = 1. The /i,f3 are real parameters and 
the rii are the Kac labels related to the negative of the highest root do = — Yli=i n i&i- The 
Virasoro central charge of the conformal Toda field theories and the conformal dimension 
of the perturbing operator V = finQe 13 airip have been computed in [21] 



„ 4£h(h + l) . , 2h , s 

where we use the effective coupling 1 < B = 20 2 j(fi 2 + 4vr) < 2. Since 2h > 1 - B/2 
is always true we are in the regime A < and expect the above mentioned complications 
with regard to renormalization to arise. To establish that, we recall first M, pL |3| 



sirB stt(2 — B) , six ... . .lai 
Vij {0) = -2 sm ^ sin 1 - >- / sin — Xl {s)x 3 (s)e~ s \ e \ , (6.3) 

see 

and deduce thereafter from (2/7) and ( 2.1 2]) 



irm 2 sm(ir/h) 
sin(7rS/2/i)sin[7r(2 - B) /2h] 



= ,.,. n rnZ."n rrwoTi • ( 6 - 4 ) 



Clearly, as < B < 2 we have (T^) > 0. Up to an overall mass re-scaling of m — > 2m, this 
agrees precisely with the results in ||12]| , which were obtained by matching the high-energy 
behaviour of the scattering matrix with a Feynman diagramatic analysis. 

It is very interesting to note that by an analytic continuation from real to purely 
complex coupling we can also reach the regime for A > for ( |6.4[ ) and observe similar 
phenomena as for the HSG-models 2 . For h = 2 this means we continue from sinh-Gordon 



1 Confusion arises sometimes due to different conventions. For instance we can relate our notations to the 
ones used in J3lj by a simple rescaling of the fields tp — ip F /\/4^r compensated by a scaling of the coupling 
constant /3 — bp/V^r- In addition, one takes the effective coupling constant to be B — 2Bf- 

2 We are grateful to Al. B. Zamolodchikov for pointing this out to us. 
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to sine-Gordon. Following for this case the argumentation of Destri and De Vega [jy], we 
relate the sinh-Gordon coupling (3 to the sine-Gordon coupling (3 via (3 — > i/3/y2 according 
to the standard conventions. Also we replace the breather mass scale m with the soliton 
mass scale m via vn? — ► 4rh 2 sin 2 ttB/2 such that we end up with the simple formula 

(T^)=7rm 2 tan|(^--) , (6.5) 

~ 2 

where A = j3 /8ir is the conformal dimension of the perturbing cos-term in the model. This 
agrees also with [ j32| . Note that in the previous argument we considered sinh-Gordon as a 
perturbed Liouville theory, whereas now we perturb the free theory rather than complex 
Liouville. Analysing (|6.5|) in more detail one observes 



Irpfl \ 



< 


for 


2n-2 
2n-l 


< A < 


♦ CO 


for 




A = 


> 


for 


2n-l 
2n 


< A < 


= 


for 




A = 



2n~l 

2n 
2n~l 

t (6-6) 



2n+l 
2n 



2n+l 

with n £ N. Note that in particular for n = 1 we have as for the homogeneous sine-Gordon 
model at A = 1/2 a transition point at which the sign changes by passing through a 
singularity. Moreover, precisely this value corresponds to the free Fermion point, which in 
this case is a very explicit example for the free Fermion picture advocated above. However, 
in this case the structure is more complicated as first of all we have an infinite number of 
such points rather than just one as in the HSG-models. In addition (T^) is not always 
divergent at these points, but can also vanish. 

6.2 Non-simply laced 

It is known, that the above mentioned complication of mass renormalization is reconciled 
if one views ATFT's in terms of dual pairs of Lie algebras. Since simply laced Lie algebras 
are self-dual, this picture does not yield anything new for that case. The dual pairs of 
non-simply laced Lie algebras are (G!p,D^), (F 4 (1) ,^ 2) ), {Bf\ A^) and {cf\ D^). 



Each algebra of these pairs allows for a description of the form ( |6.1D related to each other 
by the strong- weak duality transformation (3 — > 4-7r//3, where the untwisted algebras relate 
to the weak coupling limit. The vacuum energies associated to all non-simply laced affine 
Toda theories were stated in [^]. As in there no details were presented on how they were 
obtained, it will be instructive to show that the procedure outlined in section 2, also holds 
in this case. 

Let us first of all see what we have to expect with regard to the arguments outlined 
above and compute the Virasoro central charge and the dimension of the perturbing oper- 
ator. According to [18| we have 

c = £+12Q 2 with Q = f3p+^f, (6.7) 

with (p v ) p being the (dual) Weyl vector of the untwisted Lie algebra given by half the sum 
of the positive (co) roots. Note that when evaluating ( |6,7j ) for the simply laced case yields 
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precisely (|6.2j ), but for the non-simply laced case it differs from the expressions in |^] by 
the use of p v rather than always p. The confomal dimension of a spinless primary field 
V s (x) = e (<3+ s )£(z) i n the underlying CFT is A (a) = (Q 2 - a 2 )/2, such that the perturbing 
field / u "o^( ( g So _ ( j-j (x) has conformal dimension 

A(Pa -Q)=Pa .Q-^, (6.8) 

<3o defined as in the previous section, that is being the negative of the highest root. It will 
turn out that these dimension will always be smaller zero. We will compute the precise 
values for some concrete examples below. 

Unlike to the previous cases the expansion for the kernel (p does not appear in the 
literature, we therefore start here with the scattering matrix, which can be cast into the 



universal form [33, 34] 



Sij(±9 > 0) = exp 



T 8 J j S mh^ h t)smh(t^ H t) [K- l {t)].. e ±i; " ;r 



(6.9) 



where t?^ = (2 — B)/2h, $h = B/2H with h being the Coxeter number of the untwisted 
algebra and H its dual Coxeter number /i v multiplied by the twist of the second algebra. 
The effective coupling is now generalized to B = 2H/3 2 /(H(3 2 + Airh). The U are the 
symmetrizers of the incidence matrix of the untwisted algebra ijljj = Ulji, with t{ = 2/a 2 
and the length of long roots normalized to a 2 = 2. Also needed in fl6.9|) is the inverse 
of the (/-deformed Cartan matrix Ky(t) = (qq tj + q~ 1 q~ t i} — = with deformation 
parameters q = exp(fri?/i), q = exp(ii?#) and [hj]^ = (q Iij — q~ Itj )/(q — <Z 1 )- 

In order to evaluate the expansion for (p, we can employ once again the residue the- 
orem for a contour along the real axis closing up in the positive half of the complex 
plane encircling all poles on the imaginary axis in the upper half plane. Recalling that 
det K(t) = rise£ 4 cosh [(t + ins)/2h] cosh [(t — iirs)/2h], we know the positions of all poles 
and it follows from the integral representation ( |6.9| ) that the TBA kernels admit a series 
expansion of the form 

tp ..(0) = 16iV Res [smh(^t)sinh(L-i?jyt)^(t) i ,/detiC(t)e i * / 7r l . (6.10) 
see 

We do not have a closed formula for the cof actors K, but for the sake of our argument it 
will be sufficient here to present some examples. 

6.2.1 (G^ , )- ATFT 



Let us first compute (|6.7|) and (pjg). We carry out the computations in terms of the 



quantities of the untwisted algebra for which we have two simple roots a\ and 5.2 
normalised as a 2 = 2 = 3a\. Furthermore, the Weyl vector, its dual and the negative of 
the highest root are given by 

p = 5c?i + 3a2 , p v = 55i + 6l2 and «o = — 3c?i — 2«2 • (6-H) 
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These realizations allow to compute the quantities needed in fl6.7|) and (| 
14, 3/> v .p v = 26 and 3p.p v = 3p y .p = 8. It follows therefore 



2 + 32 



13 + 3B(B - 3) 



and A 



3B + 2 
B-2' 



that is 3p = 
(6.12) 



B(2 - B) 

Clearly for < B < 2 we have — oo < A < -1. 

To proceed further we need the (generalized) Coxeter number for this theory, which 
are h = 6 and H = 12. The symmetrizers are t\ = 3 and t2 = 1. Evaluating ( |6. 10 ) and 
reading off the first order coefficient we obtain 



Btt 



COSf (1 

where we normalized the masses to 



. (2-B)w 

r sm- — j7T— sin -5- m n mh 
8^3 ^ V o, 6 = 1,2, 



4 



mi = m cos — (1 H 

6 4 



and 



m2 = m . 



We deduce then with ( |2.12|) 



vrm 2 cosf (1 - f ) 



(6.13) 



(6.14) 



(6.15) 



12 8 

Agreement with the results in [18] is achieved by changing to the conventions used in there. 
For this one needs to re-define the effective coupling to B — ► B' = 3B /(4+B) and introduce 
a "floating" Coxeter number H' = (1 - B')h + B'h y . 



6.2.2 (F 4 (1) ,^ 2) )-ATFT 

(1) 2 

For F> ' we normalize the four simple roots to ct^ 



«2 



2a\ 



2a\ 



2. The Weyl 



vector, its dual and the negative of the highest root are in this case given by 

p = I6ai + 30a 2 + 42a 3 + 22a 4 , p v = p + 22a 4 , 
and Sq = —2ol\ — 3d?2 — 4a3 — 2on± 

such that p 2 = 39, p v -p v = 402 and p.p v = p v .p = 55. With this we find 



4 + 12 



1608 + B(331B - 1388) 



and A 



16+_B 
B-2 ' 



(2 - B)B 

Therefore — oo < A < — 8. 

For this theory we have h = 12, H = 18, t± = = 1 and £3 = t\ 
between the masses of the four particles in the theory are 

-*=2sin-(l + — ), -3 =l + 2cos- 1--), -i = 2coB— (1 
mi 4 18 mi 6 6 mi 12 

We choose the normalization such that mi = m and obtain from ( B.lOj ) 



(6.16) 



(6.17) 



2. The ratios 



6 ' 



. (2-BW ■ Btt 

£ = -8V3 |i- "~ ^ a, 6 = 1, 2, 3, 4. 



cosf (1 



_B ■ 

18 - 



m^ 



(6.18) 



(6.19) 
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Therefore with (gig) we get894/4 = ^ 



K 1 fJ-/ 



7rm 2 cosf(l-g) 
4^3 sin sin B * 



(6.20) 



is 



We can match with the formulae in by B -> B' = 45/(6 + B), H = 3(4 - B') 
mi — > m' l5 m2 — > m 3 , 7713 — ► 777,4 and 7774 — ► m 2 . 



6.2.3 (^ 1} ,^ 2) )-ATFT 



Let us now present the simplest example of the family (b\ , -A^-i) - g enera l) we choose 
for the algebra £?) the normalizations d-j = 2 for i = 1, ...,£— 1 and = 1. Then 
have 

2p = 3ai+4a2, 2/9 V = 3«i + 802 and <3o = — qi — 2a2, (6.21 
from which we compute 12/? 2 = 30, p v -p v = p 2 + 72 and /?.p v = p 2 + 4. Therefore 



we 



c= 2 + 8 
Hence —00 < A < — 2 



447 + 245(473 - 17) 
B(2-B) 



and A 



5 + 4 
B-2 



For this theory we have h = 4, H = 6, t\ = 1 and £2 = 2. The masses satisfy 



mi 
m 2 



2sm-(l + - 



and we choose m\ = m. Evaluating ( |6.10 ) we obtain now 



(6.22) 



■ (2-BW ■ Stt 



sin f(l + f) 



m z 



a, 6 



1,2, 



(6.23) 



and therefore with (2.12) 



^m 2 S inf(l + f) 

(2-B)tt „;„ Btt 
6 



(6.24) 



Defining once more B — > 5' = 45/(6+5) and 5 = 4—5' we find agreement with [jDjj l. The 
previous results also hold for the (Cg , -D3 )-theory by exchanging the roles of particles 1 
and 2, since the Dynkin diagrams of 5^ and C^p are identical up to the exchange of the 
short and the long root. 

These examples are sufficient to support the validity of the approach outlined in section 

2. 



7. Conclusions 

We used the thermodynamic Bethe ansatz to compute vacuum energies (T^) for various 
types of perturbed conformal field theories. Despite the fact, that the models considered 
exhibit different general behaviours, the assumption i)-iii), needed for the validity of the 
approximations in the TBA, hold in all cases. 
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The general behaviour of (T^) is shown to be sensitive to IR and UV-conterterms, 
whose presence can be characterized by the conformal scaling dimension A of the perturbing 
operator. In the regime < A < 1/2, realized by minimal ATFT and g^-HSG models for 
k > /i v , the quantity (T^) can be identified with the IR-counterterms needed to compensate 
the divergencies in the perturbative series expansion (jl.2j), when viewed on a cylinder. 
In contrast, in the regime 1/2 < A < 1, realized by g^-HSG models for k < h v , the 
quantity (T^) can be associated to the sum of the aforementioned IR counterterms and UV 
counterterms needed to guarantee the finiteness of the individual integrals in the expansion. 
In the models studied here these additional counterterms, when passing from A < 1/2 to 
A > 1/2 show up in a change of sign in (T^). It would be extremely interesting to verify 
this assertion by some explicit perturbative computations for the HSG-models. For the 
regime A < 0, realized here by the ATFT (simply laced as well as non-simply laced) (T^) 
constitutes a mixture of several types of counterterms, less obvious to disentangle. The 
divergence of {T^} at A = 1/2 can be attributed to the occurrence of free Fermions, for 
which such type of behaviour is well known from explicit analytical expressions. However, 
we were not able to identify the free Fermions in all g/jV-HSG models, which can be viewed 
as perturbed CFT's with A = 1/2. This needs further investigations. 
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